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Abstract 



[~^ I Supersymmetric gauge models with local horizontal symmetries are known 

to generate large flavor changing neutral current efl'ects induced by super- 
symmetry breaking D-terms. We show how the presence of a U{1) gauge 
symmetry solves this problem. We then construct a realistic gauge model 
with SU{2)h X U{1)h as the local horizontal symmetry and suggest that the 



^O \ U{1)h factor may be identified with the anomalous U{1) induced by string 

^r^ ' compactification. This model explains the observed hierarchies among the 

Q^ ■ quark masses and mixing angles, accommodates naturally the solar and at- 

mospheric neutrino data, and provides simultaneously a solution to the su- 
Qh, persymmetric flavor problem. The model can be excluded if the rare decay 

Qh! u — > 67 is not observed in the current round of experiments. 

j>^ : I. INTRODUCTION 

One of the fundamental puzzles of the standard model is a lack of understanding of the 
fermion mass and mixing hierarchies. A promising approach to resolve this puzzle is to use 
horizontal symmetries, either global and local, that transform fermions of one generation into 
another |T|]. The fact that in the limit of vanishing Yukawa couplings, the standard model 
has an enormous [SU{3)1 x 5f/(3)^ x ^f/(3)^ x 5f/(3)i x 5f/(3)^ x U{1)b-l] horizontal 
symmetry makes this approach quite plausible. 

In this paper we shall be concerned with local realization of horizontal symmetries. Local 
symmetries have certain clear-cut advantages over their global counterparts. The most 
significant difference is perhaps the strong suspicion that global (but not local) symmetries 
are susceptible to explicit violation through quantum gravitational effects. 

The presence of supersymmetry appears to put additional constraints on models with 
local horizontal symmetries |0] . For instance, if one chooses a general simple non-Abelian 
group G as the local horizontal symmetry, in the presence of supersymmetry breaking, 
the D-terms associated with G will induce a non-negligible splitting among the slepton 
and squark masses of different generations. In turn they induce unacceptable amount of 
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flavor changing neutral current effects. Since the squark mass splittings are independent 
of the horizontal gauge couplings as long as the horizontal group is simple, there is no free 
parameter that can be used to dial down these splittings to an acceptable level. This poses 
a serious roadblock to the use of local horizontal symmetries as a way to understand the 
fermion mass hierarchies in a supersymmetric context. 

If the D-term mass splittings between squarks (or sleptons) were absent (or if they were 
under control), local horizontal symmetries can neatly address the SUSY flavor problem that 
plagues the generic softly broken supersymmetric standard model. Such a solution would be 
highly desirable especially in scenarios where supersymmetry breaking is communicated to 
the MSSM sector through supergravity. The clear advantage is that no special assumption 
need be made about the Kahler potential (or the superpotential) , apart from the requirement 
of gauge invariance. Two different approaches have been adopted in the past that evade the 
aforementioned D-teim difficulty: (i) assume the horizontal symmetry to be global 0, or 
(ii) use a discrete gauge symmetry [^. In both cases there are no associated Z)-terms. 

In this paper we propose a solution to the D-term problem that would facilitate the 
use of true gauge symmetries to address simultaneously the fermion mass problem and 
the supersymmetric flavor problem^. We show that by adjoining a local U{1) symmetry 
to the existing non-Abelian horizontal symmetry G, the mass splittings between different 
generations of squarks caused by the D-terms can be brought under control. The mass 
splittings will now depend quadratically on the ratio of the two gauge couplings, which can 
be adjusted to make the FCNC effects sufficiently small. We illustrate this mechanism by 
using the example of a horizontal SU{2)h x U{1)h model. We then construct a fully realistic 
model using this group and show how the f/(l) factor may be identified with the anomalous 
f/(l) that arises in superstring compactification. We find this model to be quite predictive 
in the neutrino as well as in the quark and lepton sectors. Specifically, it supports the large 
angle atmospheric and small angle solar neutrino oscillations. The model can be excluded 
if the rare decay fi —>■ e'j is not observed in the current round of experiments. Furthermore, 
this model can easily be grand unified into an SU{5) or 5*0(10) group. 

II. SUPPRESSION OF D TERM SPLITTINGS IN THE PRESENCE OF AN 

EXTRA U{1) 

Let us consider the gauge group of the theory to be SU{3)c x SU{2)l x f/(l)y x Gh, 
where the horizontal group Gh is chosen to be SU{2)h x U{1)h as stated before. We will 
choose the matter content of the model to be the same as that of the MSSM with the 
addition of three right-handed neutrinos (denoted by z/f ) . These z/f fields are required for 
the cancellation of triangle and global SU{2)h anomalies. An immediate consequence is 
that the left-handed neutrinos will have small masses induced by the seesaw mechanism. 
The horizontal quantum numbers are chosen under the straightforward assumption that the 
particles of the first two generation belong to doublets of the SU{2)h group whereas the 
particles of the third generation are in the singlet representation. (This assignment is same 



^For a specific gauged 5*0(3) model with symmetry breaking via fundamental fields, see Ref. 



as in Ref. [Q]). Furthermore, we impose the unifiabihty condition, which means that all 16 
members of a generation have the same horizontal charges. Thus the matter content of the 
model in the standard notation is: 

{Q„, L„, <, <, e^, vl) : 2(-l) 

{Q3, ^3, ul dl el, ul) : 1(0) (1) 

where we have exhibited the SU{2)h x U{1)h quantum numbers. Here a = 1 — 2 is the 
SU{2)h index. The Higgs sector consists of the following fields: 

{iJ„, if,} : 1(0); 0, :2(1); 0, : 2(-l); x : 1(1); X:l(-1); S, : 1(0) (z = 1 - 2) . (2) 

Here ii^, Hd are the usual MSSM doublet fields, while all the other fields are singlets of the 
standard model. Given the horizontal symmetry group G, this Higgs system is the minimal 
choice that can properly break G without breaking supersymmetry. The (0, 0) fields break 
SU{2)h completely, while the (x, x) fields break U{1)h- The Si fields are necessary to allow 
cubic terms in the superpotential, a requirement if the horizontal gauge symmetry is to be 
broken in the supersymmetric limit with only renormalizble terms. 

Let us first demonstrate how the suppression of the D-term mass splittings between 
different generation squarks (and sleptons) arises in this model. In this particular model the 
problem concerns the SU{2)h -D-terms which can potentially split the masses of the first 
two generations. Note that the U{1)h -D-term will not induce mass splittings within the 
first two generations. The most general superpotential involving the Higgs fields has the 
form: 

W = /i^00 + A00^i + W'ixx, S,) , (3) 

where we wont need the explicit form of the piece W involving {xXy Si) fields. In the 
supersymmetric limit, we have (0) = (0) = V^ and (x) = (x) = V^. We shall make 
the reasonable assumption that the scale of horizontal symmetry breaking (V^, V^) is much 
greater than the scale of soft supersymmetry breaking terms, Msusy- The requirement 
of vanishing F-terms in the supersymmetric limit implies F^ = (/x,^ + AS'i)0 = and 
Fs = Xcfxp + dW /dSi = 0. Including arbitrary soft supersymmetry breaking, the scalar 
potential involving (0, 0) fields is given by: 

v=\f^^+ xsf (10P + m + |a00 + dw'/ds\' + ^{gij, + 4^?^)(i0r - \wy 

+ m||0|2 + m||0|2 + {5<^/i^00 + A^Xcf^S + H.c.} . (4) 

The coefficients in the last line of Eq. (4) (m^, m|-, B"^, A^) are all of order M|usy- Here we 
have adopted the point of view that in the absence of a symmetry, the superpartner masses 
at the Planck scale should be arbitrary. Specifically, we do not assume universality of scalar 
masses. 

Minimizing this potential (Eq. (4)) with respect of and fields and subtracting the 
two extremization conditions, we arrive at the relation: 
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^{\ct>ct> + dW'/dS) - -{qIh + Hh) + ^0/^0 + A^^S 



^j-^4>) ■ (5) 



Noting that A# + dW'/dSi = 0{MsvsyV^), Eq. (5) implies: 

\i9lH + 4^m)(l0l' - l^n ^ 'H - ml) + 0(Msusy/1^^)M|usy ■ (6) 

The contribution to the squark (or slepton) mass sphttings from the D-term is given by: 

Am| ^ Isl^M' - m - , 'l\ , H - ml) . (7) 

Note that in the absence of the U{1)h group (i.e., if gm = 0), Am? ~ m| — m^ which is 
independent of the gauge couphng and is arbitrary (i.e., anywhere between (100 GeV)^ — 
(1000 GeV)^. Ahhough these D-terms contribute to diagonal squark masses, because they 
are not universal, once Cabibbo rotation on the quark fields are made, they do contribute 
to flavor changing processes. The most stringent constraints arise from K^ — K^ mixing and 
the rare decay /i -^ ej. The Ki — Ks mass difference sets a constraint |^ (from the most 
stringent {LL)[RR) operator) [Am?/m?]6'(7 < 1 x 10~^(mg/500 GeV), where m^ denotes the 
average squark mass. The constraint arising from // — »• 67 is similar. Clearly, if gm — >■ 0, 
the D-term splittings will grossly contradict experiments^ if the squark masses are below a 
TeV. On the other hand, in the presence of the extra gauge coupling giH-, we can control the 
FCNC processes to adequate levels. For example, if g2H/ 9ih = (1/3 — V^)) which is not at 
all unreasonable, then Am? ~ (1/40 - l/200)(m? - m|). For m^ ~ (300 - 500) GeV, this is 
comfortably consistent with experiments. Note that the soft supersymmetry breaking mass 
terms m| and m|- do not run below the horizontal scale (since the and fields have masses 
of order V^). On the other hand, the masses of the squarks do run below V^ and in this 
process receive a flavor universal contribution from the gauginos. For comparable values of 
initial (Planck scale) gaugino and squark masses, the factor m^^/m,"^ is suppressed by about 
~ 1/10 because of the running. 

While we used a speciflc example to illustrate the proposed mechanism to cure the D 
term problem, its features prevail in more general contexts. For example, we could use 
alternative superpotentials (Cf. Eq. (3)) such as W = XSi^cfxp — /i^) + W or one involving 
non-renormalizable operators. The FCNC problem in the absence of U{1)h, and its solution 
via U{1)h will be identical to the case discussed above. 

III. A REALISTIC MODEL OF FERMION MASS AND MIXING HIERARCHIES 

We will now show that the model presented in the previous section can naturally explain 
the fermion mass and mixing hierarchies. We take the viewpoint that all operators consistent 
with gauge invariance are allowed in the Lagrangian. This includes non-renormalizable 
operators, which will be suppressed by appropriate inverse powers of the Planck mass. The 
coefficients of such operators will all be assumed to be of order unity. 

Let us first consider the quark sector of the theory. It is easy to see that the superpotential 
consistent with SU{2)h x U{1)h symmetry is given by: 



^For an exception where 9c arises almost entirely from the up-quark sector, see Ref. Q. 



+ ^,^^'QaUlH^x' + ^,^^'QadlH,x' . (8) 

Defining two small parameters e^ = -^ and ex = -^j we get the following hierarchical mass 
matrix for both up and the down sectors: 



Mf = Vf 



( hi^el \ 

"'12^;,,. ''■22^0 ''■23^0 

V hi^e^ hi^ ] 



(9) 



with f = u,d. The charge lepton mass matrix has an identical form as Eq. (9), as does the 
Dirac neutrino mass matrix (identify / = £, z/ for the two cases). 

The mass matrices in Eq. (9) naturally explain the fermion mass and mixing angle 
hierarchies. To see this in detail, assume that all the hij parameters are of order one. The 
mass ratios in the down-quark sector is then given by: 

ms/rrib ~ e^, md/rris ~ e^/e^ (10) 

with similar results in the up-quark and the charged lepton sectors. If we choose e^ ~ 1/7 
and e^ ^ 1/20, all the observed masses and mixing angles can be explained, with the 
coefficients hij taking values in the range (1/2 — 2). This is a tremendous improvement over 
the standard model Yukawa couplings which span some six orders of magnitude. In our 
scheme, order one differences such as in m^/m,- and ms/rrib (the two differ by about a factor 
of 3 near the Planck scale) are attributed to order one differences in the hij couplings. The 
hierarchy mi)/ nit requires moderate to large values of tan/5 ~ 10 — 40. 

The zeros in the mass matrices of Eq. (9) are corrected only at very high order, and are 
negligible. (The (1,1) entry receives a correction of order e^e^, the (1,3) and (3,1) entries 
are corrected at order Ct^e^.) The near vanishing of the (1,1) entry, along with the (anti)- 
symmetry of the (1,2) entry leads to a successful quantitative prediction for the Cabibbo 
angle; the vanishing of the (1,3) and (3,1) entries yield a relation for Vub (and Vtd)'- 



IK„ 



^Jmjms - ^Jmjm^e'"' ; |Kfe|/|Vcb| ~ ^niu/m^ ; \Vtd\ / \Vts\ ^ \lrnd/rri, . {II) 



The last two relations could serve as future tests of the model. 

Turning now to the leptonic sector, as noted, the charged lepton and the Dirac neutrino 
mass matrices have identical forms as Eq. (9). The z/f Majorana mass matrix is obtained 
from the Lagrangian: 

+ ^e'^'Ki^lA^bX' + ^^>,^A0^0^6»^e^'^x' + ^<KM,e^'e'\' . (12) 



At the level of the standard model, a bare mass term will be allowed for the u^ Majorana 
mass. We have assumed it to arise from the VEV of a field A that breaks B — L symmetry. 
When the model is embedded into a left-right symmetric or S'O(IO) framework, or if the 
B — L symmetry of the model as it stands is gauged, the Majorana masses of the v'^ fields will 
require the VEV of such a multiplet. Apart from the motivation to unify, we follow this path 
since then i?-parity violating terms will be automatically eliminated from the Lagrangian. 
Due to the intricacy of the seesaw diagonalization, we have kept the lowest non-vanishing 
terms in all entries of the Majorana v'^ matrix, which is given by: 



M,e = (A) h^elel h2el /sse^ ■ (13) 




V /i3e</,e2 
Using the seesaw formula the light left-handed neutrino mass matrix is obtained to be: 

2 / -Fiie^/e^ Fuel/el Fiae^/e^ \ 

M^^'* = - TXV ^1^4/4 F22e^" F,,e-/ ■ (14) 

^^^ V F,sel/e^ F23e^' F33 / 

Here Fij are functions of various combinations of /i^- and fij, and are expected to be of order 
one. 

It is amusing to note that the largest entry in Eq. (14) corresponds to the mass of u^^. 
The light neutrino mass hierarchy predicted in the model is m,^^ -C m^^ <^ m^^ . If we set 
(A) = 2 X 10^^ GeV, which is the supersymmetric unification scale, 771^, ^ 7 x 10~^ eV. 
This is in the right range to explain the atmospheric neutrino data via z/^ — z/^ oscillations. 
The relevant mixing angle is given by O'^^ ~ (-F23/-P22 ~ ^23/^33)^0- This angle can be near 
maximal (45°) if, for example, F23/F22 ~ 3 — 4. This appears quite plausible, given that Fij 
are non-trivial combinations of the Yukawa couplings. {-F23/F22 = {fi3f22h'^3 + fi2f2?,h^3 — 
/i2/33/i32 + /i3/23^32)/[(^i2(/23 " /22/33)]}- As for the solar ueutriuo problem, it is explained 
by small angle v^ ~ ^t MSW oscillations. The mass of v^ is of order A x niu^ ~ 10~^ eV. 
The Vf^ — Vj. mixing angle is of order e^/e^ ^ 0.02. Both parameters neatly fit the desired 
values 0. Large angle Ve oscillations are very unlikely in this model, if it is established the 
model could be excluded. 



IV. SOLVING THE SUPERSYMMETRIC FLAVOR PROBLEM 

The model presented here has a built-in solution to the supersymmetric flavor problem. 
In fact, part of our motivation to use the horizontal symmetry was to address this problem. 
Since the horizontal symmetry is local, no explicit violation is expected from quantum 
gravitational effects. As already discussed (see Sec. II), augmenting the horizontal symmetry 
group by a U{1) factor alleviates flavor violation arising from the horizontal SU{2)h D~ 
terms. Flavor violation in the squark (and the slepton) sector will however be induced, once 
effects of horizontal symmetry breaking are included. We will now show that such violations 
are not excessive and are consistent with present FCNC constraints. 

We assume that supersymmetry breaking occurs in a hidden sector and is communicated 
to the MSSM sector by supergravity. However, we do not make any special assumption about 



the Kahler potential or the superpotentiaL In particular, we do not assume that the scalar 
masses are universal or that the supersymmetry breaking trilinear A-terms are proportional 
to the superpotential Yukawa couplings. The soft scalar masses arise from general Kahler 
potential terms. For the first two generation squarks the dominant contribution arises from 

C=fd'eQiQ^-^ (15) 

•' ^^-' Planck 

where z is a hidden sector (spurion) field with nonzero F-component that breaks super- 
symmetry. Identifying F^/Mpi^^^^^. = M^^o^y, we see that the dominant masses for the first 
two generations are universal as a consequence of the SU{2)h horizontal symmetry. Non- 
universal corrections will arise from terms such as 

C= j d%Q\<PMlQb)j^ (16) 

■' ^'^ Planck 

and a similar term with is replaced by </). Compared to the dominant contribution (Eq. 
(15), these non-universal terms are suppressed by a factor e?, ~ 2 x 10~^. Since these 
corrections contribute to diagonal entries in the squark mass matrix, any FCNC effect will 
require an additional quark mixing angle (~ Oc ~ Jf^d/ms ~ 0.2). Furthermore, as noted 
earlier, these non-universal corrections get diluted by about a factor of 1/10 in the squark 
sector through RGE effects proportional to the gaugino (mainly the gluino) mass. It is 
convenient to define a set of parameters 5^2 as the ratio of the (1,2) entry of the squark mass 
matrix to the average squark mass-squared in a basis where the quark fields are physical 
0. We then estimate {5f2)LL,RR ~ (2 x 10"^) x (0.2) x (0.1) = 4 x 10"^ This is to be 
compared with the experimental limit on this quantity, {5f2)LL,RR < 1 x 10^^ valid for an 
average squark mass of 500 GeV p. We see broad agreement with experiment. 

Analogous discussion in the first two generation slepton sector leads to a prediction 
{5[2)ll,rr ^ e^^eM ^ (2 X 10-2) X (0.07) ~ 1.4 x 10"^. Here we have taken the e - /^ 
mixing angle to be Jmf./m^ ~ 0.07, appropriate to the mass matrix of Eq. (9). Note that 
unlike the squark sector, there is no significant dilution effect due to the RGE (since sleptons 
are color neutral). This number should be compared with the constraint from the present 
experimental limit on /x — > 67, which is {5{2)ll,rr ^ (4.0 x 10^'^ — 1.8 x 10^^) corresponding 
to nil — 100 GeV and for x = m'^/m'- in the range 0.3 — 3 [|1^. Although the constraint 



is satisfied, the rate for jj. ^ e'j cannot be much below the present experimental limit. We 
estimate the rate to be at most a factor of 100 below the present limit, which will soon be 
tested. 

As for the supersymmetry breaking trilinear A terms, they arise in supergravity models 
from superpotential terms such as 

C= f d'OQ.dlHd—^ . (17) 

The resulting coefficients of the trilinear scalar terms are of order MgusY- Horizontal gauge 
invariance implies that in our model, the structure of the A-terms is identical to that of the 
superpotential in Eq. (8). However, the coefficients of the A matrix are not proportional 
to the Yukawa matrix arising from Eq. (8). This non-proportionality will lead to FCNC 



processes. We estimate the parameter {Sf2)LR ~ ^'^Avd/m'^- Choosing niq = 500 GeV and 
tan/5 = 30, we obtain {5f2)LR — 3 x 10~^, which is well below the experimental limit on this 
quantity arising from K^ — K^ mixing {{5f2)LR < 3 x 10~^). A similar estimate (perhaps 
slightly smaller value, since Jme/rn^ ~ (1/3) Jrud/ms) will hold for the leptonic {S{2)lr,rl- 
For slepton masses of 500 GeV, the constraint from yu ^- 67 is {S12)lr,rl < 2 x 10~^. We 
see that the constraint is quite tight. Allowing for unknown order one coefficients (or some 
proportionality of the A terms) we conclude that /i — *> 67 cannot be much below the present 
experimental limit. Since the coefficients of j1l(^r and flRei are approximately the same, we 
expect that both helicity muons will participate in the decay, unlike the grand unification 
effects discussed in Ref. ||ll| where the decay fii -^ 6^7 is suppressed. The rate for the 



decay r ^ /X7 is estimated to be two orders of magnitude below the present limits. 

V. COMMENTS AND CONCLUSIONS 

Before concluding a few remarks are in order. 

(i) An important check for the renormalizability of the model is that the horizontal gauge 
symmetries be anomaly free. SU{2)h is automatically free of chiral anomalies. Cancellations 
of anomalies involving the U{1)h should however be ensured. With the assignment of U{1)h 
charges as given, we find the anomaly coefficients for U{1)h x [SU{2)l]'^, U{1)h x [SU{3)c]'^, 
U{1)h X [SU{2)h]'^ and U{1)h x [U{1)y]^ are all equal to -8. (We have used the GUT 
normalization for the Y quantum number.) It is then very tempting to identify the U{1)h 
as the pseudo-anomalous U{1) that arises in superstring compactification [|I^. An attractive 
aspect of this identification is that the x ^^^ X fields can play the role of the singlet fields 
for the purpose of model building 1113 and the dominant mode of supersymmetry breaking 



could be via the D-terms of the anomalous U{1) group. In that case, one of them (which 
can be chosen to be the x field) picks up a VEV of order ^Mpianck- Thus if we scale all 
higher dimensional operators by the Planck mass, we get the desired order for the e^. It 
is then clear that we must choose (0) ~ /0\ also of the same order. This identification 
will go well with the fits obtained from the fermion masses. The x could be used to break 



supersymmetry if we kept only the XX term in the superpotential as in Ref. |]T3[. We do 
not pursue this line here. It might be mentioned that anomaly cancellation can also occur 
by introducing fields which are vector-like under the standard model, but chiral under the 

U{1)h. 

(ii) The model as it stands can easily be embedded into the grand unification groups such 
as SU{^) or 5*0(10) since the horizontal quantum numbers are common to all the quarks 
and leptons that fit into a single multiplet of the above groups. As usual, for the case of 
SU{b) unification the MSSM doublets must be embedded into 5 and 5 representations of 
SU{b). For the case of S'O(IO) grand unification, additional multiplets belonging to 126 or 
16 will be needed. 

(iii) CP violation in the model can arise via complex Yukawa couplings as in the stan- 
dard model. There are additional supersymmetric source for CP violation. An interesting 
possibility is that e'/e in the K meson system can be explained through the gluino pen- 
guin. As noted earlier, the parameter {5f2)LR,RL — 3 x 10~^ in our model. If its argument 
is of order one, it is of the right magnitude to explain e'/e ~ 2.7 x 10"'^ 0]. {5f^LR,RL-i 



which arises after Cabibbo rotation, will be of order 6 x 10~^. This will lead to a neutron 
(and electron) electric dipole moment very near the present experimental limit. It should, 
however, be noted that the horizontal symmetry by itself does not fully cure the electric 
dipole moment problem, since the B parameter and the gaugino masses should have small 
phases. An interesting scenario would be one where SUSY breaking terms arising through 
the superpotential only (the A-terms) have order one phases. 

In conclusion, we have discussed a way to avoid excessive FCNC effects induced through 
D-term mass splittings between squarks and sleptons of different generations in mod- 
els with local horizontal symmetries. We have constructed an anomaly free model with 
SU{2)h X U{1)h as the local horizontal symmetry group and shown that it can lead to a 
proper understanding of the observed hierarchies among the quark and lepton masses and 
their mixings. We have shown that the model provides a simultaneous solution to the su- 
persymmetric flavor problem. Without any additional assumption this model also leads to 
a desirable pattern of neutrino masses and mixings; it supports small angle oscillations for 
the solar and large angle oscillations for the atmospheric neutrino data. Although flavor 
violation in the model is under control, it is not unobservable. The rare decay /x — > 67 is 
predicted to be near the present experimental limit. 

The work of K.S.B is supported by funds from the Oklahoma State University. R.N.M. 
is supported by the National Science Foundation grant No. PHY-9802551. 



REFERENCES 

[1] For some early attemps see: R. N. Mohapatra, Phys. Rev. D9, 3461 (1974); 

H. Georgi, Phys. Rev. DIO, 539 (1974); 

F. Wilczek and A. Zee, Phys. Rev. Lett. 42, 421 (1979) and references therein. 
[2] Y. Kawamura, H. Murayama and M. Yamaguchi, Phys. Rev. D51, 1337 (1995). 
[3] M. Dine, A. Kagan and R. Leigh, Phys. Rev. D48, 4269 (1993); 

A Pomarol and D. Tommasini, Nucl. Phys. B466, 3 (1996); 

R. Barbieri, G. Dvah and L. Hall, Phys. Lett. B377, 76 (1996); 

T. Blazek, S. Raby and K. Tobe, |hep-ph/9903340 . 



[4] D. Kaplan and M. Schmaltz, Phys. Rev. D49, 3741 (1994); 

L. Hall and H. Murayama, Phys. Rev. Lett. 75, 3985 (1995). 
[5] K.S. Babu and S.M. Barr, Phys. Lett. B387, 87 (1996). 
[6] M. Ciuchini et. al, |hep-ph/ 9808328 . 



[7] Y. Nir and N. Seiberg, Phys. Lett. B309, 337 (1993). 

[8] L. Hall and A. Rasin, Phys. Lett. B315, 164 (1993). 

[9] For recent reviews of the the neutrino data as well as its possible interpretations, see: 

K. Zuber, Phys. Rept. 305, 295 (1998) ; 

S. Bilenkii, C. Giunti and W. Grimus, |hep-ph/9812360| . 
[10] Y. Kim, P. Ko, J. Lee and K. Lee, Phys. Rev. D59, 055018 (1999). 
[11] R. Barbieri, L. Hall and A. Strumia, Nucl. Phys. B445, 219 (1995). 
[12] M. Dine, N. Seiberg and E. Witten, Nucl. Phys. B 289, 589 (1987); 

J. J. Attick, L. J. Dixon and A. Sen, Nucl. Phys. B 292, 109 (1987). 
[13] P. Binetruy and E. Dudas, Phys. Lett. B389, 503 (1996); 

G. Dvali and A. Pomarol, Phys. Rev. Lett. 77, 3728 (1996); 

R. N. Mohapatra and A. Riotto, Phys. Rev. D 55, 1138 (1997); Phys. Rev. D 58, 1138 

(1997); 

A. Farragi and J. C. Pati, Nucl. Phys. B 526, 18 (1998); 

A. Nelson and D. Wright, Phys. Rev. D 56, 1598 (1997); 

R. Zhang, Phys. Lett. B402, 101 (1997); 

N. Arkani-Hamed, M. Dine and S. Martin, Phys. Lett. B431, 329 (1998); 

K. Kurosawa and N. Maekawa, [hep-ph/9902469| . 



[14] A. Masiero and H. Murayama, |iep-ph/9903363 



K.S. Babu, B. Dutta and R.N. Mohapatra, |hep-ph/9905464 . 



10 



